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ABSTRACT 

The Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey has detected high-mass star¬ 
forming clumps with anomalous N 2 H+/HCO+( 1 - 0 ) integrated intensity ratios that are either un¬ 
usually high (“N 2 H+ rich”) or unusually low (“N 2 H+ poor”). With 3 mm observations from the 
Australia Telescope Gompact Array (ATCA), we imaged two N 2 H+ rich clumps, G333.234-00.061 
and G345.144-00.216, and two N 2 H+ poor clumps, G351.409+00.567 and G353.229+00.672. In these 
clumps, the N 2 H+ rich anomalies arise from extreme self-absorption of the HCO^ line. G333.234- 
00.061 contains two of the most massive protostellar cores known with diameters of less than 0.1 pc, 
separated by a projected distance of only 0.12 pc. Unexpectedly, the higher mass core appears to be at 
an earlier evolutionary stage than the lower mass core, which may suggest that two different epochs 
of high-mass star formation can occur in close proximity. Through careful analysis of the ATCA 


observations and MALT90 clumps (including the 
regions), we find that N 2 H+ poor anomalies arise 
abundances of N 2 H+ due to the distinct chemistry 
Subject headings: stars: formation - stars: massi 

1. INTRODUCTION 

The formation of high-mass stars (>8 Mq) is not as 
well understood as the formation of low-mass stars (e.g., 
Zinnecker & Yorke 2007). Part of the reason is that 
high-mass stars are more disruptive to their natal envi¬ 
ronment due to large accretion rates, strong winds, and 
ionizing photons. The chemistry of star-forming regions 
mostly depends on the temperature and density, but in 
high-mass star-forming regions the effects of shocks and 
photochemistry become important as well. High-mass 
star-forming regions are rarer and more deeply embed¬ 
ded compared with their low-mass counterparts, further 
adding to the difficulty in understanding the formation 
processes which alter the chemistry of the natal environ¬ 
ment of high-mass stars. 

High-mass stars form within high-mass cores (defined 
here to have a size-scale of < 0.1 pc and mass < 100 Mq) 
that arise in high-mass clumps (defined here to be ^1 pc 
and > 200 Mq). Since high-mass stars form in regions 
with large visual extinctions, observations at wavelengths 
longer than visible light are needed to determine their 
locations and study their formation processes. Contin¬ 
uum surveys of the Galactic plane have been undertaken 
at multiple wavelengths, including GLIMPSE {Spitzer 
IRAC 3.6, 4.5, 5.8, and 8.0 /im, Benjamin et al. 2003), 
MIPSGAL {Spitzer MIPS 24 and 70 /im, Carey et al. 
2009), Hi-GAL {Herschel PACS/SPIRE 70, 160, 250, 
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G333, NGC 6334, and NGC 6357 star formation 
at clump-scales and are caused by lower relative 
of H II regions or photodissociation regions. 

- astrochemistry - ISM: clouds - ISM: molecules 

350, and 500 /im, Molinari et al. 2010), and ATLASGAL 
(APEX, 870 /im, Schuller et al. 2009). The ATLASGAL 
survey was particularly successful at identifying high- 
mass star-forming clumps at a variety of evolutionary 
stages (Foster et al. 2011). 

Based on the locations of the ATLASGAL clumps, 
the Millimetre Astronomy Legacy Team 90 GHz survey 
(MALT90, Foster et al. 2011; Jackson et al. 2013; Foster 
et al. 2013) used Mopra to observe 16 line transitions to 
explore the chemical composition and kinematic informa¬ 
tion of ^3000 of these clumps. The catalog for MALT90 
sources is given in J. Rathborne et al. (in prepara¬ 
tion), which provides velocities, integrated intensities, 
and linewidths of each line for every observed ATLAS¬ 
GAL clump. Moreover, J. Rathborne et al. (in prepara¬ 
tion) classifies each source based on its infrared emission 
(using GLIMPSE and MIPSGAL) as Quiescent, Proto¬ 
stellar, H II region. Photodissociation Region (PDR), or 
Unknown (details are in Section 2.2). 

Using the first year of MALT90 data, Hoq et al. 
(2013) explored the chemical evolution of clumps and 
found several “N 2 H+ anomalies”. These sources showed 
unexpectedly extreme integrated intensity ratios be¬ 
tween N 2 H+( 1 - 0 ) and HCO+(1-0). “N 2 H+ rich” 

sources were defined to have integrated intensity ratios 
[I(N2H+)/I(HG0+)] greater than 4 and “N 2 H+ poor” 
sources have [I(N 2 H+)/I(HGO+)] less than 0.3. 

Ghemical models for star-forming regions suggest 
that extreme values of the integrated intensity ratio 
[I(N 2 H+)/I(HG 0 +)] can occur as the star-forming re¬ 
gion evolves. During the prestellar phase, low temper¬ 
atures coupled with high densities cause GO to freeze 
onto dust grains, while nitrogen-bearing molecules re¬ 
main relatively unaffected (e.g., Gharnley 1997; Bergin 
& Langer 1997). As star formation progresses, the tem¬ 
perature rises, and GO is released from the dust grains 
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and destroys N2H+ molecules through the chemical re¬ 
action N2H+ + CO ^ HCO+ + N2 (e.g., Busquet et ah 
2011). Thus it may be expected that cold, dense sources 
are N2H+ rich, while warmer H II regions and PDRs are 
N2H+ poor. 

To investigate these intensity anomalies, we used the 
Australia Telescope Compact Array (ATCA) to observe 
two N2H+ rich sources and two N2H+ poor sources. We 
discuss the classification of MALT90 sources, the ATCA 
source selection, and the data reduction in Section 2. 
We analyze the N2H+ rich sources in Section 3 and find 
that G333.234-00.061 contains two very high-mass pro- 
tostellar cores. We discuss the two observed N2H+ poor 
sources and develop a sample of all MALT90 N2H+ poor 
sources in Section 4. Finally, we summarize the findings 
in Section 5. 

2. CLASSIFICATION, SOURCE SELECTION, 
OBSERVATIONS, AND DATA REDUCTION 

2.1. MALT90 Classification of Sources 

J. Rathborne et al. (in preparation) classified all 
observed MALT90 sources as Quiescent, Protostellar, 
H II region, PDR, or Unknown. This classifica¬ 
tion was done by eye and was based on two dif¬ 
ferent Spitzer mid-infrared RGB images: (1) 3.6, 
4.5, and 8.0 pm (GLIMPSE), and (2) 3.6, 8.0, and 
24 pm (GLIMPSE+MIPSGAL). The classification is 
based solely on the emission contained within the Mopra 
MALT90 beam (38") centered at the location of an AT¬ 
LAS GAL source. An example of the classifying scheme 
can be seen in Eigure 1 of Jackson et al. (2013). The 
ATCA and MOPRA dishes have the same diameter and 
thus the ATCA primary beam is also well described by 
the MALT90 classification. 

“Quiescent” (or prestellar) sources have no obvious 
emission at 3.6, 4.5, 8.0, and 24 pm and appear as dark 
extinction features at these mid-IR wavelengths. “Pro¬ 
tostellar” clumps show a compact region of 24 pm emis¬ 
sion usually associated with more extended regions of 
enhanced 4.5 pm emission (“green fuzzies,” Chambers 
et al. 2009). “H II regions” contain high-mass stars with 
strong UV radiation, causing emission to be bright and 
extended in all bands. They are particularly bright in 
the 8 pm (due to fluorescent excitation of polycyclic aro¬ 
matic hydrocarbons) and 24 pm bands, causing these 
sources to appear yellow in the GLIMPSE+MIPSGAL 
images. “PDR” sources are at the interface of molecu¬ 
lar and ionized gas and show extended emission in all 
bands. They are especially bright in the 8 pm band and 
consist of well-defined PDR ridges. Einally, sources that 
did not fit into these four main categories were classified 
as “Unknown”. These clumps usually had no obvious 
feature in the Spitzer images or contained a combina¬ 
tion of extinction features intermixed with patchy 8 /im 
emission. 

Hoq et al. (2013) found that for the year 1 MALT90 
data, dust temperatures derived from Herschel/Yii-GAh 
continuum data correlate with this clump classification 
scheme, i.e., the temperature rises from Quiescent to Pro¬ 
tostellar to H II/PDR. This correlation provides confi¬ 
dence for the clump evolutionary stage assignment. 

2.2. Selection and Parameters of Anomalous Sources 


To investigate the N2H+ anomalies, we chose 
two N2H+ rich sources, AGAL333.234-00.061 and 
AGAL345.144-00.216, and two N2H+ poor sources, 
AGAL351.409+00.567 and AGAL353.229-I-00.672. 
Eor simplicity in the rest of this paper, we refer 
to these sources by their “G” names (G333.234- 
00.061, G345.144-00.216, G351.409H-00.567, and 

G353.229H-00.672) rather than their ATLASGAL 
(AGAL) names. 

These sources represent some of the most extreme val¬ 
ues of the [I(N2H+)/I(HC0+)] ratio and have some of the 
brightest spectral lines of the anomalous sources. The 
two N2H+ rich sources are shown in Eigure 1, and the two 
N2H+ poor sources are shown in Eigure 2. The two N2H+ 
rich sources each contain a bright region of enhanced 4.5 
pm emission (a “green fuzzy”, Eigure 3a, 3b) and com¬ 
pact 24 pm emission (Eigure 1) indicative of a protostel¬ 
lar source. The two N2H+ poor sources are dark in all 
bands (Eigure 2, 3c, 3d). Thus, contrary to prediction, 
the N2H^ rich sources contain compact heated regions, 
which should produce HCO+ and destroy N2H+, while 
the N2H+ poor sources are dark, which could indicate a 
cold source that should be rich in N2H+. Spectral en¬ 
ergy distribution fits to the Hersehel Hi-GAL data sug¬ 
gest that these clumps have masses larger than 200 Mq 
(Y. Gontreras, in preparation) and therefore are likely 
to produce at least one high-mass star >8 Mq (Jackson 
et al. 2013). 

Kinematic distances based on the MALT90 molecu¬ 
lar velocities for G333.234-00.061 and G345.144-00.216 
are 4.76lo'29 1-8610'57 kpc respectively (S. Whitaker 

et ah, in preparation). The other two dark clumps, 
G351.409H-00.567 and G353.229+00.672, are poorly de¬ 
fined by kinematic distance because they lie at longitudes 
close to the Galactic center. However, G351.409H-00.567 
and G353.229H-00.672 are coincident, both spatially and 
in velocity, with the star-forming complexes NGC 6334 
and NGG 6357, respectively, which are reported in Rus- 
seil et al. (2012) to have distances of 1.7 ± 0.3 kpc and 
1.9 ± 0.4 kpc. In this paper, we adopt these distances. 

Eollowing Hoq et al. (2013), the dust temperature 
and mass column density for each clump can be esti¬ 
mated from the Hersehel/Hi-GAL Survey (A. Guzman 
et al. in preparation). These parameters were de¬ 
rived by fitting a gray body model with a single 
temperature to the Hersehel Hi-GAL 160, 250, 350, 
and 500 pm maps. Each map was smoothed to the 
500 pm resolution, allowing for the dust tempera¬ 
ture estimates to be applicable on the scale of 36". 
The central pixel values of the temperature maps for 
G333.234-00.061, G345.144-00.216, G351.409-H00.567, 
and G353.229-H00.672 are 21, 18, 22, and 30 K respec¬ 
tively, and the central pixel for the mass column density 
is 0.75, 0.17, 0.074, and 0.075 g cm“^, respectively (A. 
Guzman, private communication). Using the distances 
for each source as discussed above, this suggests masses 
within the 36" Hersehel beam of approximately 3100, 
no, 35, and 45 M© for G333.234-00.061, G345.144- 
00.216, G351.409-H00.567, and G353.229+00.672 respec¬ 
tively. The clump for each source is extended and much 
larger than this size (as seen by the ATLASGAL emission 
in Eigure 3), and thus is much more massive; these mass 
estimates simply give the approximate masses within the 
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Fig. 1.— Spitzer three-color image (3.6 pm. (blue), 8.0 pm (green), and 24 /rm(red)) of the two N 2 H+ rich sources overlaid with MALT90 
HCO+(1-0), N2H+(1-0), HCN(l-O), and HNC(l-O) contours. Left: G333.234-00.061 with contours shown for [3, 6, 9, 15, and 24] x cr, 
where a = 0.36 K kms“^. Right: G345.144-00.216 with contours shown for [3, 5, 7, 9, and 13] x cr, where a = 0.37 K kms“^. 
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Fig. 2.— Spitzer three-color images (3.6 pm (blue), 8.0 pm (green), and 24 pm (red)) of the two N 2 H+ poor sources overlaid with 
Mopra MALT90 HGO+(1-0), N2H+(l-0), HGN(l-O), and HNG(l-O) contours. Left: G351.409+00.567 with contours shown for [3, 6, 
12, 18, and 24] x cr, where cr = 0.49 K kms“^. For clarity, some contours are numbered based on their signal-to-noise, and thicker 
contours indicate more statistically significant contours. Note that HGO+ and HGN show significantly stronger emission than N 2 H+. 
Right: G353.229+00.672 with contours shown for [3, 6, 9, 12, 18, and 24] x cr, where cr = 0.33 K kms“^. The green artifacts on the 
bottom left of this image are due to saturation in the MIPS 24 pm band. 


ATCA 3 mm primary beam. Moreover, Herschel wave¬ 
lengths are primarily sensitive to the cold component of 
gas, and thus a hot compact core (as probed by ATCA) 
that is much less massive than the cold clump mass can 
be completely hidden in the Herschel bands. 

For the two N2H+ rich sources, G333.234-00.061 and 
G345.144-00.216, Spitzer IRAC identifies protostars 
that are likely the hottest sources within the dark clump, 
suggesting the temperatures of 21 and 18 K, respectively, 
are lower limits to the core temperatures. The other two 
sources, G351.409+00.567 and G353.229+00.672, have 
no obvious protostellar sources, and thus the cores in 
these clumps may be much colder since dense environ¬ 


ments can shield from external radiation. 

2.3. Observations and Data Reduction 

Using the ATGA we observed four molecular lines, 
HCO+(1-0), HNG(l-O), 13CS(2-1), and N2H+(l-0), and 
two 2-GHz wide continuum bands (at 90 and 93 GHz) 
over the course of four days: from 24 August to 27 Au¬ 
gust 2013. The array configuration was H168, which pro¬ 
vided an angular resolution of ^2.2" at these frequen¬ 
cies. The Gompact Array Broadband Backend correla¬ 
tor setup for the observations was GFB 64M-32k, which 
provides 64 MHz zoom bands (i.e., a velocity range of 
^200 kms“^) with 32 khz (0.11 kms“^) spectral resolu- 
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Fig. 3.— Spitzer three-color image (3.6 pm (blue), 4.8 pm (green), and 8.0 pm (red)) of the four ATCA sources observed in this paper 
overlaid with contours from APEX ATLASGAL 850 pm continuum observations, a) G333.234-00.061, with contours shown for [3, 10, 
30, 70] X (7, where cr = 54 mjybeam“^; b) G345.144-00.216, with contours shown for [3, 6 , 10, 14] x cr, where cr = 58 mjybeam“^; c) 
G351.409+00.567, with contours shown for [3, 5, 7, 9] x cr, where a = 72 mjybeam”^; and d) G353.229+00.672, with contours shown for 
[3, 6 , 10, 14] X cr, where cr = 61 mJybeam”^. 


tion along with two bands of 2 GHz bandwidth of contin¬ 
uum. For data reduction, we used MIRIAD (Sault et ah 
1995) and reduced the data in the standard manner. All 
four sources were observed every day, and the data from 
each day were combined. The phase calibrator used for 
G333.234-00.061 and G345.144-00.216 was 1600-44 and 
for G351.409+00.567 and G353.229+00.672 was 1742- 
289. For bandpass calibration, 1921-293 was used for 3 
of the 4 days, while 1253-055 was used for the other day. 
Uranus was used as a flux calibrator, and the flux cal¬ 
ibration is expected to be accurate within ^20%. Rest 
frequencies for HGO+(1-0), HNG(l-O), and ^^GS(2-1) 
were adopted from the Splatalogue database^ and were 
taken to be 89.188526, 90.663572, and 92.49427 GHz re¬ 
spectively. A more accurate N2H+(l-0) rest frequency 
of 93.173772 GHz was adopted from Daniel et al. (2006). 
The continuum observations were combined, resulting in 
a center frequency of 91.5 GHz. 

The shortest baseline for this configuration is 61 m, 
which, at 91.5 GHz, makes these observations insensi¬ 
tive to extended emission larger than about 11". Since 
these observations are only sensitive to structure smaller 
than 11" and the sources have extended emission for the 
brightest molecular lines (i.e., HGO+, HNG, and N2H+), 
very large negative sidelobes are evident in our images. 
Trying to suppress such large sidelobes proved to be fu¬ 
tile. Instead, we generated “dirty” images using natu¬ 
ral weighting and identified the brightest emission peaks. 
These peak intensity regions were then “cleaned” using 
MIRIAD’s clean algorithm. 

Gontinuum emission was detected toward the two pro- 
tostellar sources, G333.234-00.061 and G345.144-00.216; 
this emission was subtracted from the molecular line ob¬ 
servations. All ATGA fluxes reported in this paper have 
been corrected for primary beam attenuation. 

Also discussed in this paper are 22.2 GHz continuum 
observations of G333.234-00.061 taken with the ATGA 
on 5 May 2014 in the 1.5D array configuration. These 
data were reduced in a similar manner as the 91.5 GHz 
observations, using 1253-055 as the bandpass calibrator, 
1646-50 as the phase calibrator, and 1934-638 as the flux 
calibrator. Fluxes at 22.2 GHz are accurate to within 

® Spectral line data were taken from the Spectral Line 
Atlas of Interstellar Molecules (SLAIM) (Available at 
http://www.splatalogue.net). (F. J. Lovas, private commu¬ 
nication, Remijan et al. 2007). 


TABLE 1 

MALT90 Spectral Lines 


Species 

Main transition 

Frequency (GHz) 

N 2 H+ 

J = 1- 0 

93.173772 


J = 2-1 

92.494303 

H 

41a 

92.034475 

GH 3 GN 

1 

II 

91.985313 

HG 3 N 

J = 10 - 9 

90.979020 

13C34S 

J = 2-1 

90.926036 

HNC 

J = 1- 0 

90.663572 

HG^^CGN 

J=10-9, F = 9- 8 

90.593059 

HCO+ 

J = 1- 0 

89.188526 

HGN 

J =1-0 

88.631847 

HNGO 

JKa,Ki, = 4o,4 - 3o,3 

88.239027 

HNGO 

JKa,Ki, = 4i,3 - 3i,2 

87.925238 

G 2 H 

jv = 1 - 0 , J = l-i 

F = 2 - 1 

87.316925 

HNi^C 

J = 1-0 

87.090859 

SiO 

J = 2-1 

86.847010 

h13co+ 

J = 1-0 

86.754330 


Note. — Bold rows indicate spectral lines that 
were also observed with ATCA for the four sources 
in this paper. 


-10%. 

Additionally for these four sources, we present the 
spectra from the MALT90 survey (J. Rathborne et al. 
in preparation). These lines are summarized in Table 1. 

3. N 2 H+ RIGH SOURGES 
3.1. G333.234-00.061 
3.1.1. Continuum 

The continuum map for G333.234-00.061 is shown in 
Figure 4 and shows two very bright 3.3 mm continuum 
cores in the field. The southern core (henceforth MMl) is 
significantly brighter than the northern core (henceforth 
MM2). A weaker third continuum core may be present 
between the two main continuum peaks. Neither MMl 
or MM2 is resolved. The continuum observations have a 
synthesized beam size of 274 (calculated as the geometric 
mean of the synthesized beam’s full width at half maxi¬ 
mum (FWHM) major and minor axes), corresponding to 
an upper limit for the core sizes of —0.06 pc. 

We used MIRIAD task maxf it to locate the intensity 
maxima; the results place MMl at R.A. (J2000.0) = 
16^19“5H29 and deck = -50°15'1474 and MM2 at 
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R.A. (J2000.0) = 16^19“50!94 and decl. = -50n5'10'/5. 
These peaks are separated by 5'/l or a projected distance 
of 0.12 ± 0.012 pc (where the error only conies from the 
uncertainty in distance). Since these cores have similar 
velocities (as discussed in Section 3.1.2) and are associ¬ 
ated with the same MALT90 clump, their separation is 
probably not much larger than this projected distance. 
The clump velocity is a better indicator of the kinematic 
distances than the actual core velocities because MMl 
and MM2 may have peculiar velocities within the clump. 
If we were to assume the projected distance is the same 
as the radial distance, the separation of the cores would 
increase to 0.17 pc. However, since we are unable to de¬ 
termine the radial distance with the given data, we will 
simply refer to the separation as a “projected distance of 
0.12 pc” in the rest of the paper. 

MM 2 is associated with Spitzer IRAC emission in all 
four bands (only 8.0 pm is shown in Figure 4) and with 
Spitzer MIPS at 24 pm. MMl is not detected in any 
of the IRAC bands. The resolution for MIPS at 24 pm 
(^ 6 ") is insufficient to resolve the two cores (separated 
by 571), though the 24 pm emission is centered on MM2 
and is not elongated toward MMl. This 24 /im compact 
emission and presence of a green fuzzy (Figures 1 and 
Figure 3) indicates that MM2 is likely protostellar and 
MMl is colder and at a younger age; however, MMl is 
likely protostellar as well since it is coincident with a class 
II (radiatively pumped) methanol maser, as discussed in 
Section 3.1.3. 



51.8s 51.6s 51.4s 51.2s 51.0s 50.8s 16hl9m50.6s 


RA (J2000) 

Fig. 4.— The ATCA 3.3 mm continuum image for G333.234- 
00.061, corrected for the primary beam. White contours and 
color-scale are shown for the 3.3 mm continuum, with contours 
shown for [-3, 3, 5, 7, 10, 14, 20] X ( 73 .3mm, where ( 73 . 3 mm = 
0.17 mJybeam“^. Blue contours are the Spitzer IRAC 8.0 pm 
with contour levels [200, 300, 500, 1000] x 1 MJysr“^. The white, 
black, and green x marks show locations of class I methanol masers 
for 36, 44 (Voronkov et al. 2014), and 95.1 GHz (Ellingsen 2005), 
respectively. The red, cyan, and yellow x marks indicate 6.7 GHz 
class II methanol (Caswell et al. 2011), 22 GHz water (Breen et al. 
2010), and 1.7 GHz OH masers (Caswell 1998), respectively. 


Since MMl and MM 2 are unresolved, we approximate 
the integrated flux densities to be the same as the peak 
flux densities integrated about the synthesized beam, giv¬ 
ing integrated flux densities of 19.0 and 15.1 mJy for 
MMl and MM2 respectively. The total flux density of 


the continuum is 47.5 mJy (calculated by summing the 
flux in a polygon around the emission), giving a leftover 
integrated flux density (which is dominated by the con¬ 
necting dust lane) of 13.4 mJy. 

To calculate the dust masses of the cores, we use the 
standard equation from Hildebrand (1983) 


Md= 


( 1 ) 


where is the observed integrated flux density, D is the 
distance to the source, is the dust opacity, and Bj^(Td) 
is the Planck function given at dust temperature To. As 
discussed in Section 2 . 2 , Herschel data show that the 
cold clump mass in the ATCA beam is large, so the de¬ 
rived dust temperature at the clump-scale (36" resolu¬ 
tion) is 21 K and is likely a lower limit oi To for these 
protostellar sources. The 3 mm continuum emission is 
not typically dominated by the clump or the central pro¬ 
tostar but rather by the thermal dust emission from the 
surrounding dust cocoon at a radius of ^1000 to 10000 
AU. Dust temperatures in a typical high-mass hot molec¬ 
ular core at these radii vary from ^30 K to 120 K (Os¬ 
orio et al. 1999). For this paper, we adopt To = 100 K. 
This temperature is likely high for MMl since the lack 
of IRAC emission indicates it may be much colder than 
MM 2 . Assuming that oc and /3, the dust emis- 
sivity index, is 1.5, we scaled 6 : 1 . 3 mm = 0-90 cm^ g“^ 
(Ossenkopf & Henning 1994, assuming thin ice mantles 
and gas density of 10 ^ cm“^) to 6 : 3.3 mm = 0.22 cm^ g“^. 
We assume a gas-to-dust mass ratio (GDR) of 100 to es¬ 
timate the core mass. From these values, we calculate 
gas masses of 36 Mq and 29 Mq for MMl and MM2 
respectively. The remaining material of the complex at 
this temperature has a mass of 26 Mq, suggesting the 
total mass of the core complex to be ^100 Mq. 

Estimating errors of these masses is difficult since the 
parameters are highly uncertain. If our assumed values 
are the most probable, we can estimate errors for our 
sources. We assume a 20% error in F 3 3 and 6 : 1.3 mm- As 
discussed in Section 2 . 2 , the kinematic distance with un¬ 
certainties is D = 4.76lo.29* errors of the other 

parameters, we assume they are uniformly distributed 
in the following ranges: GDR G (70,150), P G (1,2), 
and T G (20 K, 120 K). Our mass estimates with errors 
for MMl and MM2 are thus 361^5 Mq and 29^.12 ATq. 
However, our parameters were selected because they give 
conservatively low mass estimates. For example, if we 
change the Td, one of the most dependent variables for 
estimating the mass, from 100 K to 50 K (with the same 
range for each parameter), the masses would be 74^40 

and 59^30 Mq respectively. For the rest of this section 
we will use the warmer masses of 36 Mq and 29 Mq for 
MMl and MM2. Given that the parameters for Equation 
1 were chosen to give conservatively low masses, these 
masses are likely lower limits. The ratio between the 
masses of MMl and MM2, however, is much less uncer¬ 
tain since all parameters in Equation 1 are likely simi¬ 
lar for both sources except FI and possibly T^. There¬ 
fore, MMl is at least 25% more massive than MM2 since 
MMl’s flux is 25% higher and it is likely at a colder tem¬ 
perature. 

In the infrared dark cloud SDC335.579-0.272, Peretto 
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et al. (2013) found the highest 3 mm luminosity ever 
seen for a protostellar core {F^ 2 mm = 101 at a dis¬ 
tance of 3.25 kpc) for a core diameter <0.1 pc. Assuming 
Td = 50 K and / 1 : 3.2 mm = 0.087 cm^ g“^ (using f3 = 2), 
Peretto et al. estimated a core mass of 545 Mq. Indeed, 
if we assume the same dust opacity and temperatures for 
MMl and MM2, both of these cores would have masses 
around 200 Mq. Peretto et al. compared masses for 
about 30 protostellar cores, standardizing their masses 
with the same dust opacity law; such a comparison sug¬ 
gests that MMl and MM2 are perhaps the second and 
third most massive protostellar cores known for diame¬ 
ters <0.1 pc. Regardless of the exact values of k and Td, 
MMl and MM2 are two of the most massive protostellar 
cores known, and they are forming in close proximity to 
each other. In this paper, we have assumed a higher Td 
and lower [3 compared to most of the studies mentioned 
in Peretto et ah, which lowers our mass estimates. 

For these cores we calculated the beam-averaged H 2 
column density, ^"( 112 ), using 

pp 

iV(H2) = „ y -(2) 

where is the peak flux density, rrip is the mass of 
a proton, /i is the mean molecular weight per hydrogen 
molecule, and is the solid angle subtended by the 
beam. Given of 19.0 and 15.1 mJybeam”^ for MMl 
and MM2, respectively, and assuming /i = 2.8 (Kauff- 
mann et al. 2008) and a gas-to-dust mass ratio of 100, 
we determine a A^(H 2 ) of 4.6 x 10^^ and 3.7 x 10^^ cm“^ 
for MMl and MM2. 

We calculated the mean volume density of MMl and 
MM2 using the expression 


n(H2) 


Mg 

lirripV^ 


( 3 ) 


given the gas mass M^, mass of a proton m^, and an 
assumed spherical volume V = |7rr^, where r is the ra¬ 
dius of the core. Since the cores are unresolved, we use 
a radius of half the synthesized beam (i.e., a radius of 
0.03 pc), which provides a lower limit of the mean vol¬ 
ume density. The n(H 2 ) for MMl and MM2 is 4.6 x 10^ 
and 3.7 X 10^ cm“^ respectively. 

Thus far, we assumed that all of the 3.3 mm continuum 
flux arises from thermal dust emission and have neglected 
any contamination from free-free emission. Additional 
ATCA observations at 22.2 GHz (1.35 cm) detected both 
MMl and MM2. These cores are resolved with a synthe¬ 
sized beam of 1733 x 1707 (~0.03 pc resolution) and have 
peak fluxes of 0.30 and 0.39 mJybeam”^ for MMl and 
MM2 respectively. If we smooth these observations to 
the same resolution as the 91.5 GHz continuum obser¬ 
vations, the peak fluxes are 0.66 mJybeam”^ and 0.70 
mJybeam”^ for MMl and MM2 respectively. Assuming 
that flux scales as oc z/^, the spectral index a is 2.4 
and 2.2 for MMl and MM2 respectively. Although op¬ 
tically thick free-free emission can have spectral indices 
up to two, typically hypercompact (HC) H II regions 
(sizes less than ^0.03 pc. Franco et al. 2000) have a ^ 1 
(e.g., Cyganowski et al. 2011) which may be due to non- 
uniform gas density (Hoare et al. 2007). Dust emission 
from high-mass protostellar objects has a spectral index 


(a > 2, with typical values of 3 to 4.5. Toward MMl and 
MM2, thermal dust emission is certainly significant since 
a is larger than 2, but a is lower than expected for gray 
body dust emission at a single dust temperature. 

If we make the conservative assumption that all the 
91.5 GHz continuum emission is from free-free emission, 
the brightness temperature of the free-free emission can 
be calculated via 


n = 


2v‘^kVtb 




( 4 ) 


where is the solid angle subtended by the synthe¬ 
sized beam (3702 x 1793), 0 is the beam filling fac¬ 
tor, Tg is the electron temperature, and is the optical 
depth. Although MMl and MM2 could be optically thin 
at 91.5 GHz, a is 2.4 and 2.2 for MMl and MM2 re¬ 
spectively. If free-free emission indeed dominates the 3.3 
mm continuum flux, these spectral indices can only be 
explained by optically thick free-free emission. For op¬ 
tically thick free-free emission, the beam filling factor is 
therefore (j) = T\)jT^ — where Qg is the solid angle 

subtended by the source. Although Tg = 10^ K is typi¬ 
cally assumed in literature, electron temperatures in the 
Galaxy have been reported to be as low as ^5000 K (e.g., 
Quireza et al. 2006), and we have used this conservative 
value for our calculation to maximize 0. The measured 
T 5 from the peak flux densities at 91.5 GHz are 0.48 and 
0.38 K for MMl and MM2 respectively. For optically- 
thick free-free emission, this results in a maximum 0 of 
9.5x10“^ and 7.6x10“^ for MMl and MM2, suggesting 
that the sources have diameters of 07024 and 07021 re¬ 
spectively. This indicates a very small maximum size of 
^100 AU for optically thick free-free emission. An H II 
region this small likely only occurs for a brief period of 
a massive protostar’s life - this is extremely unlikely to 
be seen for two sources. Thus, we conclude that free-free 
dominant emission is unlikely and that most of our emis¬ 
sion at 91.5 GHz comes from dust. The same analysis 
of the 22.2 GHz observations also finds a maximum size 
of ^100 AU for an optically thick free-free source that 
dominates the observed emission. If we assume all the 
emission at 22.2 GHz is optically thin free-free emission 
with a fiat spectral index {a = 0 ), it would contribute 
less than 5% to the total flux at 91.5 GHz. 

A mix of thermal dust emission and free-free emission 
is still possible. Observations in other star-forming re¬ 
gions have shown that although 91.5 GHz is dominated 
by dust, frequencies lower than ^40 GHz can be domi¬ 
nated by a thermal jet (Brooks et al. 2007) or free-free 
emission from a (HC) H II region (Zhang et al. 2014), 
though a HC H II region is not likely, as explained in 
the previous paragraph. Brooks et al. (2007) showed 
that when using only 25 and 88 GHz continuum obser¬ 
vations, the spectral index of the young stellar object 
IRAS 16547-4247 is a = 1.6. However, between 88 GHz 
and 250 GHz, a = 4.7. Thus, the fact that MMl and 
MM2 have spectral indices less than expected could sim¬ 
ply be due to a mix between free-free emission from a 
thermal jet and dust emission. Since we only have two 
continuum points at high angular resolution, we cannot 
estimate the free-free contamination of the dust emission. 
However, since the 3.3 mm emission appears to be dom¬ 
inated by dust, we neglect other possible components. 
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3.1.2. Molecular lines 

Figure 5 shows the 16 lines observed in the MALT90 
survey (resolution of 38") toward G333.234-00.061, 
which provide a comparison of the clump (~1 pc 
size scale) molecular properties to that of the cores 
(<0.1 pc size scale) probed by the ATCA. With the 
Mopra 22 m single dish, HCO+(1-0) and ^^CS( 2 - 1 ) 
are undetected, HNC(l-O) is marginally detected, and 
N 2 H+( 1 - 0 ) is strongly detected. The N 2 H+ line is 
broad, indicating a large velocity dispersion or multi¬ 
ple velocities components. Many other lines are de¬ 
tected or marginally detected: C 2 H(l- 0 ), H^^CO+( 1 - 0 ), 
HC 3 N( 1 - 0 ), HN13C(1-0), HNCO( 4 o, 4 - 3 o, 4 ), and SiO( 2 - 
1 ). The fact that HCO+ is undetected and HNC is 
marginally detected but and HN^^C are de¬ 

tected indicates that the optically thicker HCO+ and 
HNC lines are likely self-absorbed. 

Integrated intensity maps of the molecular lines ob¬ 
served with the ATCA are shown in Figure 6, with strong 
detections for all four of the observed molecular lines 
(HCO+, HNC, i^CS, and N 2 H+). For the 3 strongest 
lines, HCO+, HNC, and N 2 H+, the presence of extended 
emission has caused large sidelobes. Contours have been 
set in Figure 6 to show minimal sidelobes, which causes 
the integrated intensity maps not to show some real emis¬ 
sion (e.g., MMl has definite detections for all lines). 
Thus, Figure 6 only shows the locations of the strongest 
emission. 

HCO+ and HNC are especially strong for MM 2 while 
the N 2 H+ is primarily concentrated in an elongated 
structure south of both cores. ^^CS, the most optically 
thin of the probed high-density tracers, traces well the 
two continuum cores and is also especially strong toward 
MM2. 

The fact that ^^CS is much brighter for the lower mass 
core MM2 is intriguing. Estimating the optical depth 
for ^^CS without observations of other CS isotopologues 
is difficult, especially since the compact emission is not 
resolved. Therefore, we take an analytical approach in 
order to discern why ^^CS is brighter for MM2. The 
critical density of ^^CS( 2 - 1 ) depends on temperature and 
is a few times 10 ^ cm“^, which is well below the lower 
limit of the densities of the cores (few times 10 ^ cm“^. 
Section 3.1.1). Therefore the lines are thermalized for 
both cores, causing the excitation temperature, Tex, to 
be approximately equal to the gas kinetic temperature. 
The ^^CS brightness temperature can be characterized 
via T 5 = 0Tex(l — In the optically thin limit, Ti) = 

0TexTjy, and for the optically thick limit T 5 = 0Tex- Since 
^^CS is brightest for MM2, the ^^CS emission for MMl, 
when compared to MM2, (1) has a lower optical depth 
(in the optically thin limit), ( 2 ) has a lower excitation 
temperature, (3) has more self-absorption, and/or (4) 
comes from a smaller region (i.e., has a smaller (p). Since 
column densities as calculated by dust (Section 3.1.1) are 
higher than MMl, (1) would suggest that in the optically 
thin limit, ^^CS in MM2 must be more abundant unless 
the MM2 excitation temperature is much higher. (2) 
and (3) would both suggest that MMl is colder than 
MM2. (4) suggests that MMl is more compact. Thus, 
this analysis suggests that MMl could be colder and/or 
more compact. 


Figure 6 (right) shows selected spectra for the molec¬ 
ular lines HCO+, ^^CS, and N 2 H+ for MMl (SI, at the 
continuum peak), MM 2 (S 2 , at a location to best dis¬ 
play the outflow wings), and the average spectrum of a 
small area near S 2 . The N 2 H+ spectra are contaminated 
by negative sidelobes (this is particularly evident in SI), 
causing the fluxes to be inaccurate. 

The N 2 H+ spectrum at SI (Figure 6 , bottom right) 
suggests two velocity components, though careful analy¬ 
sis of the spectral cube indicates that a negative sidelobe 
lies between these two “components” which may have er¬ 
roneously created the two peaks. Since ^^CS is optically 
thinner than N 2 H+, ^^CS is more likely to trace the core 
while N 2 H+ is more likely to trace the clump. The fact 
that ^^CS is slightly redder than the N 2 H+ line may in¬ 
dicate that the core is moving relative to the clump. A 
single Caussian fit to the ^^CS line gives an amplitude 
A = 45 ± 5 mJybeam”^ (or 1.1 ±0.14 K), a central ve¬ 
locity V = —84.7 ± 0.27 kms“^, and velocity linewdith 
Av = 4.6 ± 0.63 kms“^. 

The HCO+ spectrum is extremely self-absorbed in 
both locations. Emission at the systemic velocity is ab¬ 
sent and only redshifted emission is detected. The HCO+ 
spectrum (as well as HNC, not shown) also shows absorp¬ 
tion against the continuum at the locations of the self¬ 
absorption. Although it is difficult to measure the ab¬ 
sorption against the continuum due to the sensitivity of 
the observations, the depth of the continuum absorption 
is comparable to the peak of MMl’s continuum emission 
of 18.9 mJybeam”^. Thus, the HCO+ line is optically 
thick. 

The spectra for S 2 (Eigure 6 , top right) show two broad 
HCO+ outflow wings with a deep self-absorption fea¬ 
ture splitting the red and blue wings. N 2 H+ and ^^CS 
emission peaks are between the HCO+ outflow wings at 
similar velocities of the HCO+ self-absorption feature. 
Similar to SI, the S 2 velocity offset between ^^CS and 
N 2 H+ may indicate that the core is moving relative to the 
clump. Moreover, there is also HCO+ absorption against 
the continuum between the outflow wings, approximately 
at MM2’s continuum peak (13.6 mJybeam”^). The ^^CS 
and N 2 H+ spectra show marginal evidence for two veloc¬ 
ity components. A ^^CS spectrum averaged over the cir¬ 
cle just southeast of MM2 is shown in the middle right 
panel of Eigure 6 . The two velocity components were 
fitted with two Gaussians. The velocity of the brighter 
component is v = — 86.0 kms“^, while that of the fainter 
component is v = —89.8 kms“^. Since the optical depth 
should be larger at positions approaching MM2, and the 
two distinct Gaussian components become more obvious 
at positions displaced from MM2, two velocity compo¬ 
nents seem to be a more probable explanation of the 
double ^^CS spectral peaks rather than self-absorption 
in this transition. 

Although the lack of HCO+ emission at the systemic 
velocity of MMl could conceivably be explained not by 
self-absorption but rather by the interferometer resolving 
out large scale emission, we tested this idea and found 
that it is unlikely. If we smooth the ATCA observations 
to the size of the Mopra beam, the resulting HCO+ line 
profile (which is not Gaussian due to the self-absorption 
feature) has a peak flux of 0.07 K and a velocity of - 
82.3 kms“^. This brightness temperature is roughly the 
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Fig. 5.— MALT90 (resolution of 38'') observations of G333.234-00.061 for 16 observed lines (J. Rathborne et al., in preparation). The 
spectra have been smoothed using a Savitzky-Golay kernel. The consensus velocity, VpsR: is indicated at the top of the figure and is shown 
as a reference line in each panel. Each panel shows different observed line transitions with the red curve showing the fit to the baseline and 
the green curve showing the fit to the line. The line fit is only shown if the signal-to-noise is larger than three, and the corresponding fit 
parameters are shown in the top-left of the corresponding panel. is the peak line intensity in K, V is the center velocity and AV is 

the linewidth in kms“^, and S/N is the signal-to-noise of the line fit. 


same as the rms noise for the smoothed Mopra spectrum 
in Figure 5. If the lack of HCO+ emission was indeed 
a case of resolving out structure, the peak brightness 
of the missing flux would be much higher than 0.07 K 
and would be spatially distributed at scales larger than 
11 " (the largest structure the ATCA observations are 
sensitive to); thus, the single dish brightness would be 
much higher than 0.07 K and would be easily detected 
by Mopra. Moreover, the fact that there is HCO+ ab¬ 
sorption against the continuum and that N 2 H+ does not 
show the typical 1:5:3 intensity ratios between the three 
main hyperfine peaks shows that the central region is 
optically thick, conducive to self-absorption. Also, as 
already discussed, the single-dish detections of isotopo- 
logues (e.g., HN^^C, H^^CO^) at intensity ratios much 
higher than the abundance ratios demonstrate 

very large optical depths. Therefore, these observations 
in conjunction with the MALT90 observations indicate 
that the HCO+ is deeply self-absorbed. 

The S 2 self-absorption feature occurs over a velocity 
space of 14 kms“^ (-97 to -83 kms“^. Figure 6 , top 
right), and the flux of the ATCA observations is approx¬ 
imately zero in the entire primary beam within this ve¬ 
locity range (not shown). This is an exceptionally broad 
for an absorption feature that is spatially distributed 
throughout the clump. In order to create such an ab¬ 
sorption feature, there must be optically thick cold gas 


over a velocity range of 14 kms“^. Possible mechanisms 
to produce such a broad absorption line include infall of 
cold gas on the clump scale or extreme saturation of a 
very thick line. 

Figure 7 shows the red and blue HCO+ components 
of G333.234-00.061. The redshifted and blueshifted 
velocity wings of MM2 spatially overlap, indicating a 
bipolar outflow closely oriented along the line of sight. 
A filament-like structure connects this outflow feature 
through MMl to an HCO+ source toward the southeast. 

3.1.3. Masers 

Masers are seen throughout G333.234-00.061, as seen 
in Figure 4. Here, we briefly discuss some background 
about these maser types, and then we infer the evolu¬ 
tionary state and some properties of MMl and MM2. 

Four masers are commonly associated with high- 
mass star-forming regions: two types of methanol 
(GH 3 OH) masers, water masers, and hydroxyl (OH) 
masers. Methanol masers are categorized into two dif¬ 
ferent classes. Glass I methanol masers are collisionally 
pumped (e.g., Leurini et al. 2004) while class H methanol 
masers are radiatively pumped (Gragg et al. 1992). Glass 
I masers are often offset from the nearest continuum peak 
(e.g., Kurtz et al. 2004) and found in shocked regions, 
particularly those associated with molecular outflows, 
(e.g., Voronkov et al. 2006; Gyganowski et al. 2009; Gan 
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Fig. 6 .— G333.234-00.061 ATCA integrated intensity maps of the molecular lines HCO+(1-0), HNC(l-O), ^^CS(2-1), and N 2 H+( 1 - 
0) overlaid on a grayscale map of the 3.3 mm continuum. Reference contours are shown for [-4, -3, 3, 4, 5, 6, 7, 8, 9, 12, 15] x / for 
/hco+ ~ 0.075, /hnc = 0.05, /iscg = 0.022, and = 0.22 Jybeam“^ kms“^. The open black beam (top panel, top right) represents 

the continuum resolution while the filled blue beam (bottom left) represents the panel’s molecular line resolution. Lines were integrated 
based on the following velocity ranges: HCO+: -113.2 to -72.0 kms“^, HNC: -98.5 to -73.3 kms“^, ^^CS: -95.3 to -79.8 kms“^, and 
N 2 H+: -99.3 to -76.2 kms“^. Selected spectra are shown on the right for the positions shown in the N 2 H+ map, with the red X marks 
showing spectra at the corresponding pixel and the red circle showing the area used to present an averaged spectrum. The spectra 

have been Gaussian smoothed and show a reference dashed line at a velocity of -90.5 kms“^. The central ^^CS spectrum is fit with 
two Gaussians (indicated in blue, with the sum in red). The redshifted (right) Gaussian has an amplitude A = 87 ± 7 mjybeam“^ (or 
2.2 ± 0.18 K), central velocity v = —86.0 ± 0.12 kms“^, and linewidth Av = 2.8 ± 0.29 kms“^ and the blueshifted (left) Gaussian has 
A = 46 ± 9 mJybeam“^ (or 1.1 ± 0.23 K), v = —89.8 ± 0.18 kms“^, and Av = 1.7 ± 0.43 kms“^. 
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Fig. 7.— The ATCA HCO+(1-0) integrated intensity ob¬ 
servations showing the red and blue components of G333.234- 
00.061. MM2 has a strong outfiow almost along the line of 
sight. Contours are shown for [-4, -3, 3, 4, 5, 6, 7, 8, 9, 12, 
15] X 0.065 Jybeam ^ kms The open black beam (top right) 
represents the continuum resolution while the filled blue beam (bot¬ 
tom left) represents the resolution for the HCO+ line. 


et al. 2013). Although typically found within high-mass 
star-forming regions, class I masers have also been re¬ 
ported toward supernova remnants (e.g., Pihlstrom et al. 


2014). Class II methanol masers have thus far only 
been found associated with high-mass star-forming re¬ 
gions (e.g., Xu et al. 2008). Class I masers are suggested 
to be associated with star-forming regions at an earlier 
evolutionary stage (e.g., infrared dark clouds) than class 
II masers (e.g., Ellingsen 2006). Class I masers also can 
occur after class II masers due to, e.g., expanding H II re¬ 
gions (Voronkov et al. 2010), but MMl and MM2 appear 
to pre-date these evolutionary stages. 

Like class I methanol masers, water masers are also 
collisionally pumped and associated with shocks and out¬ 
flows (e.g., Elitzur et al. 1989), and like class II methanol 
masers, OH masers are radiatively pumped (Cragg et al. 
2002). Water masers are suggested to appear before 
OH masers (Eorster & Caswell 1989), but after methanol 
masers (Breen et al. 2007). 

In G333.234-00.061 all four masers (methanol class I 
and class H, water, and OH) are detected. Their loca¬ 
tions are indicated in Eigure 4. Both cores are associ¬ 
ated with a class H methanol maser (6.7 GHz, Caswell 
et al. 2011), while only MMl contains a class I methanol 
maser centered on the source (95.1 GHz, Ellingsen 2005). 
Additional class I methanol masers at 36 and 44 GHz 
(Voronkov et al. 2014) are scattered throughout the re¬ 
gion rather than being centered on a particular source. 
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which may indicate active outflows. Only MM2 is asso¬ 
ciated with water (22 GHz, Breen et ah 2010) and hy¬ 
droxyl (1.7 GHz, Gaswell 1998) masers. Although many 
studies (e.g., Forster & Gaswell 1989; Ellingsen 2006; 
Breen et al. 2007) suggest that the appearance of masers 
from earlier to later evolutionary stages (often with over¬ 
lap) follows the sequence “Class I Methanols Class H 
Methanols Waters OH”, exceptions are frequently seen 
(e.g., Voronkov et al. 2006; Cyganowski et al. 2012). Such 
inconsistencies are likely due to the complexity of high- 
mass star formation (e.g., large clusters of stars with a va¬ 
riety of masses contained within diverse environments). 
Nevertheless, if masers generally follow this evolutionary 
trend, MM2 is likely more evolved due to its lack of a 
class I methanol maser and MMl’s lack of water and OH 
masers. 

Although class I masers are typically associated with 
outflows, MM2’s strong outflow does not appear to be 
associated with a class I maser. The outflow may have 
broken through most of the circumstellar material and 
is no longer producing collisions capable of exciting the 
class I maser. MMl, on the other hand, could have an 
outflow that has yet to break through its natal environ¬ 
ment, and thus produces shocked gas capable of exciting 
a class I methanol maser. 

As discussed in Section 3.1.2, from ^^GS(2-1) we derive 
a core velocity of —84.7 kms“^ for MMl, and two core 
velocities of —86.0 kms“^ and —89.8 kms“^ (with the 
latter the brightest) for MM2. For MMl, the 95.1 GHz 
class I methanol maser has a velocity of —87.4 kms“^ 
(Ellingsen 2005) and the class H methanol maser has a 
velocity of-85.3 kms“^ (Caswell et al. 2011). Eor MM2, 
the class H methanol, water, and OH masers have veloci¬ 
ties of-91.9 kms“^, -88 kms“^, and -84 kms“^ respec¬ 
tively (Caswell et al. 2011; Breen et al. 2010; Caswell 
1998). As expected, the radiatively pumped (class H 
methanol) maser for MMl is closer to the velocity of the 
core than the collisionally pumped (class I methanol) 
maser, strengthening the hypothesis that the class I 
methanol maser for MMl is due to an outflow. MM2 
is more complex, perhaps due to its two velocity com¬ 
ponents. The class H methanol and water masers are 
closer to the higher velocity component, while the OH 
maser is closer to the lower velocity component. The 
difference between the velocities of the two radiatively 
pumped (class H and OH) masers may indicate multiple 
cores in MM2. 

3.1.4. Two high-mass cores at different evolutionary stages 

Erom the preceding discussion, MM2 is likely at a later 
evolutionary stage than MMl for the following reasons: 

(1) MM2 has Spitzer IRAC emission while MMl does 
not, suggesting much warmer temperatures for MM2. 

(2) MM2 has a large outflow, while MMl does not have 
a definitive outflow. 

(3) ^^CS has stronger emission on MM2 even though 
MMl is brighter in the continuum, suggesting MMl 
could be colder. 

(4) The detected maser species for MMl and MM2 indi¬ 
cate that MM2 is likely at a later evolutionary stage. 

These cores are most likely protostellar because they 
seem to lack free-free emission (see Section 3.1.1), have 
masers that are expected to be absent in prestellar cores. 


and the fact (for MM2) that the compact mid-IR colors 
indicate a deeply embedded source. The ratio of the in¬ 
tegrated flux densities of MMl and MM2 is 1.25; since 
MMl is likely colder, it should be at least 25% more mas¬ 
sive than MM2. The projected separation between MMl 
and MM2 is only 0.12 pc; therefore two of the most mas¬ 
sive protostellar cores known are forming in close prox¬ 
imity at two different evolutionary stages, with the more 
massive core at an earlier evolutionary stage. 

During the high-mass star formation process, more 
massive cores are thought to both collapse and accrete 
faster and thus evolve faster (e.g., Zinnecker & Yorke 
2007). Therefore, it is surprising that two high-mass 
cores can form within close proximity with the younger 
core more massive than the older core. We offer the fol¬ 
lowing possible explanations for these observations: 

(1) MM2 may have started to form at a much earlier time 
than MMl. This would suggest that multiple epochs of 
massive star formation can happen in close proximity. 

(2) In the past, MM2 may have had a higher accretion 
rate than MMl, and thus evolved more quickly. How¬ 
ever, MMFs accretion rate may now be higher, but it 
has yet to become as evolved as MM2. 

(3) The relationship between free-fall time, tjj, and den¬ 
sity, p, is t// oc and this timescale likely indicates 

the approximate collapse time for forming cores MMl 
and MM2. MM2’s initial environment may have had 
a higher density and thus collapsed and evolved quicker. 
However, MMl’s initial gas reservoir may have been more 
massive and has recently been accreted to make the core 
more massive than MM2. 

(4) The apparently more massive core MMl might con¬ 
sist of two or more unresolved lower mass cores. Its ear¬ 
lier evolutionary state would then be easily explained 
since lower mass cores should evolve more slowly. Indeed, 
for intermediate mass cores in the GGD27 complex, with 
higher angular resolution observations, Eernandez-Lopez 
et al. (2011) resolved an apparently more massive core 
into two lower mass cores separated by only 1700 AU, 
well below our resolution of ^10000 AU. 

(5) MM2 could actually be more massive. MMl could 
have a higher temperature than MM2 causing underes¬ 
timates of the gas mass, but we believe this as unlikely 
given the discussion in this section. Moreover, the cores 
could have different dust opacities due to different 
chemistry or densities. 

(6) Eor core accretion models, radiative and mechanical 
feedback mechanisms cause the star’s final mass to be less 
than the prestellar core (e.g.. Tan et al. 2014). Such feed¬ 
back may have been acting longer on MM2 than MMl, 
causing MM2 to be less massive than MMl even though 
MM2 was initially a more massive prestellar core. The 
possibility is unlikely since the evolutionary difference 
between the two cores is not drastically different (both 
are likely protostellar). 

3.2. G345.144-00.216 
3.2.1. Continuum 

The 3.3 mm continuum map for G345.144-00.216 (Eig- 
ure 8) reveals three unresolved cores (<273 or 0.02 pc). 
The brightest core and the dimmest core (henceforth 
MMl and MM2 respectively) are both coincident with 
Spitzer IRAC and MIPS 24 pm emission (8.0 pm is 
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shown), with the brightest IRAC emission centered on 
MM 2 . The southwest core, MM3, is devoid of any IRAC 
or 24 fim MIPS emission. The peak fluxes for MMl, 
MM 2 , and MM3 are 1.5, 0.80, and 0.94 mJybeam”^, 
respectively, and these fluxes were used to calculate the 
masses. We assume the same values for and the 

gas-to-dust mass ratio ( 0.22 cm^ g“^ and 100 respec¬ 
tively) and reproduce the same error analysis as Section 
3.1.1. For cores with temperatures of 100 K, we find 
the masses to be O.dStoJo^ 0-23lo.i5: and 0.28lo.i8 ^0 
for MMl, MM 2 , and MM3 respectively. Assuming a 
colder temperature of 50 K, the core masses approxi¬ 
mately double to 0.921 q 73 0.48lo!38 5 0 -^ 6 ^o .45 ^0 

for MMl, MM 2 , and MM3 respectively. MM3’s mass 
may be much higher since it lacks IRAC emission, sug¬ 
gesting it is colder than the other two cores. MMl and 
MM 2 are also seen with MIPS GAL at 24 fim while MM3 
is undetected. 

The MIRIAD task maxfit was used to locate the 
peak intensities, placing MMl at R.A. (J2000.0) = 
17^05“36.447!29 and deck = -41°22'03'.'8, MM 2 at R.A. 
(J2000.0) = 17^05“36!78 and deck = -41°22'07'.'6, and 
MM3 at R.A. (J2000.0) = 17^05^36^08 and deck = 
—41°22'1378. No known masers are found in the clump. 



37.5s 37.0s 36.5s 17h05m36.0s 


RA (J2000) 

Fig. 8 .— The ATCA 3.3 mm continuum image for G345.114- 
00.216, corrected for the primary beam. White contours and color- 
scale are shown for the 3.3 mm continuum, with contours shown 
for [-3, 3, 4, 5, 6, 8] x crs.smm, where crs.smm = 0.17 mJybeam“^ 

. Blue contours are the Spitzer IRAC 8.0 p,m with contour levels 
of [200, 300, 500, 1000, 1800] x 1 MJysr“^. There are no known 
masers in the clump. 


3.2.2. Molecular lines 

Figure 9 shows the 16 lines from MALT90, with the 
N 2 H+ line the brightest. The HNC spectrum has two 
obvious peaks that may represent two velocity compo¬ 
nents or a self-absorption feature. The HCO+ spectrum 
also hints at two peaks with about the same velocities as 
HNC. The (1-0), C 2 H (1-0), (1-0), HC 3 N 

(1-0), HCN (1-0), and HN^^C (1-0) lines appear to be 
detected or marginally detected. 

The G345.144-00.216 ATCA integrated intensity maps 
are shown in Figure 10 . HCO+, HNC, and N 2 H+ were 
detected while ^^CS was undetected. The three com¬ 


pact cores unexpectedly lack strong line emission. Self¬ 
absorption may suppress the ground state ( 1 - 0 ) line in¬ 
tensity below detectable levels; such deep self-absorption 
has been seen toward other high-mass clumps, e.g., in 
the HCO+ and HNC ( 1 - 0 ) lines toward the infrared 
dark cloud G028.23-00.19 (Sanhueza et ak 2013). The 
self-absorption in G345.144-00.216 could be similar to 
what we see in G333.234-00.061, except that G333.234- 
00.061 has a strong outflow delineated by optically thin¬ 
ner HCO+ emission. 

To explore the possibility of self-absorption, we ana¬ 
lyzed spectra throughout the clump. On the right of Fig¬ 
ure 10 , two example spectra are shown: one at the N 2 H+ 
emission peak (SNP, Spectra at N 2 H+ Peak) and one 
near MM 2 (SC 2 , Spectra at Core 2 ). Toward the SNP 
position, the lines are detected at this position at the 
same velocity of approximately -15 kms“^. The spectra 
for SC 2 show that near MM 2 , there are two apparent ve¬ 
locity components in the N 2 H+ spectrum whose average 
velocity is approximately -15 kms“^. The HCO+ and 
HNC line profiles, on the other hand, are offset in ve¬ 
locity from the N 2 H+ velocity peak and are also slightly 
offset from each other. If self-absorption is not occurring, 
there must be four distinct velocity components at this 
location (two indicated by N 2 H+ and one each indicated 
by HCO^ and HNC). However, it is extremely unlikely 
that these three molecules each trace a different velocity 
component, especially since they all trace one component 
at about -15 kms“^ at the SNP N 2 H+ peak. Therefore, 
these lines are likely self-absorbed with a central velocity 
at approximately -15 kms“^. Only the blue components 
for HCO^ and HNC survive, with HCO^ appearing to be 
slightly more self-absorbed. Moreover, spectra centered 
on the intensity peak of all three cores (Figure 11) show 
almost no HCO^ or HNC emission (with MM 2 having 
the only detection at ^18 kms“^) and weak N 2 H+ emis¬ 
sion for MM 2 and MM3. These results are consistent 
with Sanhueza et ak ( 2012 ), which found that clump op¬ 
tical depths are typically slightly higher for HCO+ than 
HNC, while the optical depths for N 2 H+ are over an or¬ 
der of magnitude smaller. 

G333.234-00.061’s self-absorption was much more ob¬ 
vious since its strong outflow shows broad velocity wings 
in the spectra that were not self-absorbed. No outflow 
was detected in G345.144-00.216, resulting in the non¬ 
detection of most of the observed lines toward the dust 
cores. These results lead to an interesting conclusion: if 
no lines are detected toward a core, it does not necessar¬ 
ily mean the molecule is absent. Therefore, observations 
of less abundant (e.g., isotopologues) or more highly ex¬ 
cited transitions are needed to accurately trace these ex¬ 
tremely high column density environments. 

4. N 2 H+ POOR SOURCES 
4.1. G351409+00.567 and G353.229+00.672 

The ATCA detected no 3.3 mm continuum sources to¬ 
ward G351.409+00.567 and G353.229+00.672, and the 
literature reports no masers. The Icr 3.3 mm con¬ 
tinuum sensitivity for these two sources are 0.12 and 
0.16 mJybeam”^ respectively. If we assume undetected 
cores are no more than four times brighter than the rms 
noise within a given beam size, /^ 3 . 3 mm = 0.22 cm^ g“^, 
a gas-to-dust mass ratio of 100 , and a cold core dust 
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Fig. 9.— MALT90 (resolution of 38'') observations and fit results of G345.144-00.216 for 16 observed lines (J. Rathborne et al., in 
preparation). Caption is the same as Figure 5. Two line fits are shown for HCO+(1-0) and HNC(l-O), but only the fit parameters of the 
strongest line is shown. 


temperature of Td = 10 K, we conclude that it is likely 
that there are no cores smaller than the beam more mas¬ 
sive than 1.4 Mq and 1.9 Mq in G351.4094-00.567 and 
G353.2294-00.672 respectively. However, it is possible 
there are cores that are more extended than the synthe¬ 
sized beam but have extended emission below the sensi¬ 
tivity limit. 

Molecular lines from MALT90 are shown for 
G351.409-f00.567 and G353.229-f00.672 in Fig¬ 
ures 12 and 13 respectively. G351.4094-00.567 and 
G353.2294-00.672 have integrated intensity ratios 
[I(N2H+)/I(HG0+)] of 0.07 and 0.09 respectively. Each 
source has strong G 2 H emission which is commonly 
considered a PDR tracer (e.g., Jansen et al. 1995). 

Figures 14 and 15 show integrated intensity maps of 
the three lines (HGO+, HNG, and N 2 H+) detected to¬ 
ward both clumps. In G353.229-f00.672, the HNG, 
N 2 H+, and HCO+ line emission is spatially coincident, 
but in G351.409-f00.567, each line arises from a different 
location. The lack of spatial coincidence of the lines for 
G351.4094-00.567 could be artifacts due to the interfer¬ 
ometer being insensitive to large-scale structure. 

To measure the integrated fluxes of the HGO+ and 
N 2 H+ lines throughout the moment maps, we masked 
the negative pixels and summed the positive pixels based 
on the red-dashed regions outlined in Figures 14 and 15. 
Since line fluxes are not coincident for G351.4094-00.567, 
the integrated fluxes were drawn from different areas of 
the map where the line fluxes are concentrated. For 
G351.409+00.567 we find integrated flux densities of 


0.31 Jy kms“^ and 0.65 Jy kms“^ for HCO+ and 
N 2 H+ respectively. For G353.229+00.672 we find in¬ 
tegrated flux densities of 3.7 and 2.0 Jy kms“^ for 
HCO+ and N 2 H+ respectively. On these small angular 
scales, the [I(N 2 H+)/I(HCO+)] ratios are about 2 and 
0.5 for G351.409+00.567 and G353.229+00.672 respec¬ 
tively, which is much less anomalous than the large-scale 
MALT90 line observations (0.04 and 0.12 respectively). 
If these line fluxes are diluted to the MALT90 beam 
(38"), the fluxes would all be less than 0.4 K kms“^, 
and thus all lines would either be marginally detected 
or undetected. Since MALT90 reports significant detec¬ 
tions for HGO+ (and HNG) for G351.409+00.567 and 
G353.229+00.672 (8.3 ± 0.3 and 4.8 ± 0.2 K kms-^ re¬ 
spectively, J. Rathborne et al. in preparation), the ma¬ 
jority of the HGO+ line emission of these sources and 
the presence of the N 2 H+ poor anomaly occurs at large 
scales rather than small scales. 

4.2. MALT90 and N 2 H^ Poor Sources 

Since HGO+ is usually much more optically thick 
than N 2 H+, self-absorption cannot explain N 2 H+ poor 
sources, as it does the N 2 H+ rich sources. The ATCA 
data suggest that sources only appear N 2 H+ poor on 
large angular scales. Thus, to properly investigate 
the N 2 H+ poor sources, it is best to use large-scale 
(MALT90) observations for interpretation. 

We create an N 2 H+ poor sample based on the MALT90 
catalog (J. Rathborne et ah, in preparation) with the 
following criteria: 
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37.5s 37.0s 36.5s 36.0s 17h05m35.5s 

RA (J2000) 


Fig. 10.— G345.144-00.216 ATCA integrated intensity maps of 
the molecular lines HCO+(1-0), HNC(l-O), and N2H+(l-0) on top 
of a grayscale map of the 3.3 mm continuum. Reference contours 
are shown for [-5, -4, -3, 3, 4, 5, 6, 7, 8, 9, 12, 15] x / for /hco+ ~ 
0.025, /hnc = 0.025, and = 0.06 Jybeam“^ kms“^. The 

open black beam (top panel, top right) represents the continuum 
resolution while the filled blue beam (bottom left) represents the 
panel’s molecular line resolution. Lines were integrated based on 
the following velocity ranges: HCO+: -22.1 to -11.9 kms“^, HNC: 
-22.1 to -11.9 kms“^, and N 2 H+: -26.0 to -7.1 kms“^. Selected 
spectra (Gaussian smoothed) are shown on the right for the posi¬ 
tions shown in the N 2 H+ map, with a reference dashed line drawn 
at -14.4 kms“^. 


(1) Catalog sources are only kept if they either 

have (a) no detection in N 2 H+ but a detection in 
HCO+, or (b) if they have integrated intensity ratios: 
I(N2H+)/I(HC0+)] + < 0.2, where is the error of 

the integrated intensity ratio. These cutoffs reduce our 
sample to about 20% of the total MALT90 sources. 

(2) For sources that had no detection in N 2 H+ but 
a detection in HCO+, we only consider sources with 
I(HCO+) > 2.5 K kms In order to keep sources with 
I(N 2 H+)/I(HC0+)] < 0.2,1(HCO+) < 2.5 K kms-^ con¬ 
strains I(N 2 H+) values to be less than 0.5 K kms“^. 
0.5 K kms“^ is just above the typical MALT90 inte¬ 
grated intensity detection limit, and thus this criteria 
selects sources with I(N 2 H+)/I(HCO+)] < 0.2. 

(3) In order to only retain the clumps with the most reli¬ 
able I(HCO+) measurements, sources were removed with 
an I(HCO+) signal-to-noise ratio < 6. 

(4) Sources were removed that lie within a Galactic longi¬ 
tude of — 2° < / < 4°. These sources lie within the central 



Fig. 11.— G345.144-00.216 ATCA spectra (Gaussian smoothed) 
are shown for the peak intensity pixels of MMl, MM2, and MM3, 
with a reference dashed line drawn at -14.4 kms“^. 

molecular zone of the Galaxy where clumps are subject 
to different conditions and chemical processes than those 
in the rest of Galaxy. 

(5) Sources were removed when the fine structure lines 
of N 2 H+ fall out of the frequency range of the MALT90 
bandpass. The velocity cutoff depends on the sky fre¬ 
quency; sources were typically rejected when their sys¬ 
temic velocity was less than -93 kms“^. 

(6) Sources were removed that have HCO^ linewidths 
> 15 km s“^. These sources comprised of ^10% of the re¬ 
maining sample and are therefore anomalous themselves. 
Such sources may poorly represent the N 2 H+ poor sam¬ 
ple here. 

(7) The remaining 127 sources were visually inspected, 
and three additional sources (AGAL010.356-00.149_B, 
AGAL332.191-00.047_S, and AGAL354.838+00.372_S) 
were removed due to low signal to noise that caused 
significant difficulty in providing good fits to the N 2 H+ 
spectra. 

The evolutionary classifications (as discussed in Sec¬ 
tion 2.1) for the 124 remaining N 2 H+ poor clumps con¬ 
sisted of 41 PDRs, 36 H II regions, 5 Protostellar, 8 
Quiescent, and 34 Unknown sources. Of the total 3566 
clumps MALT90 observed, there were 362, 938, 796, 689, 
and 781 sources classified as PDRs, H II regions, Pro¬ 
tostellar, Quiescent, and Unknown, respectively. The 
probability that a random MALT90 clump is an N 2 H+ 
poor source based on our criteria is 11.3%, 3.8%, 0.6%, 
1.1%, and 4.4% for PDRs, H II regions, Protostellar, 
Quiescent, and Unknown respectively. Immediately it is 
evident that these N 2 H+ poor clumps are primarily asso¬ 
ciated with regions with high ionization (i.e., PDRs and 
H II regions, with the former very evident) as well as 
the clumps with Unknown classification. The Unknown 
clumps come from diverse environments, so we cannot 
draw any general scientific conclusions about the chem¬ 
istry from these clumps. 

These 124 N 2 H+ poor clumps are especially concen¬ 
trated in two regions: the G333 complex (Figure 16, 
showing 16 sources) and NGG 6357 (Figure 17, showing 
26 sources including G353.229+00.672). Immediately it 
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Fig. 12 .— MALT90 (resolution of 38") observations and fit results of G351.409+00.567 for 16 observed lines (J. Rathborne et al., in 
preparation). Caption is the same as Figure 5. 


is evident from these Figures that all N 2 H+ poor sources 
are located on or very near H II regions. 

Why does NGC 6357 contain so many (about one fifth) 
of these N 2 H+ poor clumps? One of the reasons is due 
to the fact that it is one of the closest (^1.9 kpc) high- 
mass star-forming regions in the Galaxy. This close 
proximity allows clumps to be resolved that would be 
blended at farther distances (Figure 17 shows that many 
N 2 H+ poor clumps are spatially near each other) and 
causes the clumps to be brighter due to the inverse square 
law. On the other hand, NGG 6334 (Figure 18, show¬ 
ing ten N 2 H+ poor clumps) is also located at a close 
distance (^1.7 kpc), has a similar size, and is brighter 
than NGG 6357 but contains almost a factor of three 
fewer N 2 H+ poor clumps. The main difference between 
NGG 6334 and NGG 6357 is that NGG 6357 has a giant 
shell and the N 2 H+ poor clumps are coincident with this 
shell. The high-mass star-forming region G333 is also rel¬ 
atively close (a typical clump has a kinematic distance of 
3.2 kpc, S. Whitaker et ah, in preparation) and also has 
the majority of its N 2 H+ poor clumps bordering shells. 
These shells are typically the photodissociated interfaces 
between molecular and ionized gas. Thus, this analy¬ 
sis suggests N 2 H+ poor clumps are particularly found in 
PDR shells. 

Although N 2 H+ poor clumps are primarily associated 
with clumps classified as H II regions and PDRs, there 
remain five N 2 H+ poor clumps classified as protostellar 
and eight as quiescent. These 13 sources are all located in 
projection against bright mid-IR Spitzer emission, and 
all but three sources have at least a marginal detection 


in G 2 H (a PDR tracer). For the three sources unde¬ 
tected in G 2 H, one source had a noisy spectrum and the 
other two were relatively weak sources (I(HCO^) < 3), 
and if observations were deeper, they may have had G 2 H 
detections as well. The classification scheme is subject 
to errors due to the projection of multiple sources along 
the line of sight. Moreover, it can often fail for distant 
sources due to foreground contamination. It is entirely 
possible that these 13 sources are in fact associated with 
PDRs or H II regions within the 38" Mopra beam. 

Since N 2 H+ poor clumps in H II and PDR regions 
are associated with large amounts of radiation, these 
clumps are likely heated to higher temperatures than 
the “quiescent” or “protostellar” clumps. In higher tem¬ 
perature gas (warmer than the sublimation temperature 
of ices), GO is released from grains into the gas phase. 
Since chemical reactions with GO both create HGO+ and 
destroy N 2 H+ (as discussed in Section 1), warmer re¬ 
gions are expected to have a larger HGO+/N 2 H+ abun¬ 
dance ratio. However, at least in some sources such as 
the two observed ATCA clumps G351.409+00.567 and 
G353.229+00.672, the anomaly is no longer apparent at 
the core scale. At this scale, the high-density material 
may be at a colder gas temperature due to shielding of 
incoming radiation allowing N 2 H+ to survive. 

5. SUMMARY 

Based on MALT90 data, Hoq et al. (2013) identified 
N 2 H+ anomalous sources: sources with either high or 
low integrated intensity ratios between the N 2 H+(l- 0 ) 
and HCO+(1-0) lines. New ATGA observations of two 
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Fig. 13.— MALT90 (resolution of 38") observations and fit results of G353.229+00.672 for 16 observed lines (J. Rathborne et al., in 
preparation). Caption is the same as Figure 5. 


N 2 H+ rich clumps and N 2 H+ poor clumps show that: 

• In the two observed N 2 H+ rich sources, the N 2 H+ 
rich anomaly is due to extreme self-absorption of 
the HCO+(1-0) line. Rather than an unusually 
high N 2 H+ to HCO+ abundance ratio, the HCO+ 
emission is suppressed because it is extremely op¬ 
tically thick. 

• The ATCA observations of the two observed N 2 H+ 
poor sources suggest this emission is anomalous on 
large (clump) scales rather than small (core) scales. 
The higher density of cores may help in shielding it 
from radiation, allowing N 2 H+ to survive on small 
scales. 

• N 2 H+ poor clumps are primarily associated with 
PDRs and H II regions, especially at locations of 
shell structures. At high temperatures, CO is more 
abundant, creating HCO+ and destroying N 2 H+. 

• Given the criteria in Section 4.2, NGC 6357 
has more N 2 H+ poor star-forming clumps than 
any region surveyed by MALT90. It contains a 
plethora of these sources due to its close distance 
and the presence of a giant PDR shell. 

The N 2 H+“rich” clump G333.234-00.061 is an un¬ 
usual star-forming region. Specifically: 

• It contains two of the most massive protostellar 
cores known separated by a projected distance of 


only 0.12 pc. These cores have estimated masses of 
at least 36 Mq and 29 Mq and diameters less than 
^0.06 pc. 

• Even though class I methanol masers are con¬ 
sidered to indicate the location of outflows, the 
strongest outflow (ejected from MM2) is not asso¬ 
ciated with a class I methanol maser, while the core 
without a definite outflow, MMl, is associated with 
this maser. MMl’s class I methanol maser might 
indicate an outflow forming in MMl and generat¬ 
ing strong collisions as it starts to break through 
its natal environment. 

• The more massive protostellar core, MMl, appears 
to be at an earlier evolutionary stage than MM2 
(since MM2 has IRAG emission, a large and strong 
outflow, stronger ^^GS emission, and later-stage 
masers). For two high-mass protostars that formed 
in the same molecular clump, it is unexpected to 
have the lower mass protostar evolving faster than 
the higher mass protostar. We suggest several in¬ 
terpretations that might explain this apparent dis¬ 
crepancy. 
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Fig. 14.— G351.409+00.567 ATCA observations. Caption is the Fig. 15.— G353.229+00.672 ATCA observations. Caption is 
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the location of G333.234-00.061, an N2H+ rich source observed with the ATGA in this paper. The bright H II complex is G333. 
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